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LONG  TERM  GOALS 

The  primary  long  term  objective  of  this  project  is  to: 

•  determine  a  fast  and  accurate  inversion  method  to  estimate  bottom  properties  in  shallow 
water. 


OBJECTIVES 

The  objectives  of  this  year’s  work  (FY12)  included: 

•  continued  study  of  a  new  deterministic  low  frequency  (LF)  geoacoustic  inversion  (G.I.) 
method  recently  featuring  the  minimization  processor  (Tolstoy,  ’  10  and  ’12); 

•  demonstration  that  horizontal  arrays  can  be  successfully  used  for  G.I.  with  the  new  method 
(presented  at  the  fall  ASA); 

•  preparation  the  LF  G.I.  method  to  be  applied  to  ranges  beyond  2  km  (KRAKEN  was  inte¬ 
grated  into  the  processing  software  in  preparation  for  the  D.  Knobles  data); 

•  initiation  of  efforts  for  the  consideration  of  frequencies  above  100Hz  at  close  range  (leading 
to  consideration  of  such  frequencies  at  longer  ranges;  this  was  sugggested  by  K.  Becker  this 
past  spring).  We  have  already  seen  that  at  ranges  beyond  1km  and  at  freqs  above  100Hz 
inversion  can  be  problematic  if  all  geoacoustic,  environmental,  and  geometric  parameters 
are  unknown. 

APPROACH 

The  LF  method  performs  an  exhaustive  search  through  a  finite  parameter  space.  We  have  seen 
previously  that  some  parameters  are  inextricably  paired,  e.g.,  source  range  rge  and  water  depth 
D.  Taking  account  of  such  relationships  can  help  to  reduce  the  search  space,  particularly  for 
frequencies  below  100Hz.  We  also  note  excellent  sensitivity  for  our  low  frequencies  to  some 
geoacoustic  parameters  such  as  sediment  sound-speed  in  the  underlying  half-space  chsp  depending 
on  frequency  /  and  rge. 

For  the  minimization  approach  proposed  here  'Pmin( p)  at  the  vector  parameter  combination  p  we 
require  that  for  p  in  our  search  space  only  the  minimum  processor  value  (over  /)  be  retained.  Thus, 
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all  frequencies  considered  at  that  point  will  have  processor  values  at  least  as  great  as  the  final  value. 
Consequently,  (1)  the  frequency  most  sensitive  to  a  parameter  dominates  the  search  automatically, 
and  (2)  resolution  has  improved  with  major  sidelobe  reduction.  The  method  is  intended  primarily 
for  geoacoustic  inversion  methods  (where  the  signal-to-noise  levels  are  high). 

For  the  simulations  we  shall  continue  to  consider  one  single  sediment  layer  scenario  (each  defined 
by  a  linear  sound-speed  profile  and  constant  density)  over  a  half-space  (constant  sound-speed  and 
density)  each  at  multiple  frequencies  (25  to  100Hz,  although  recent  efforts  have  begun  investiga¬ 
tions  for  100  to  200  Hz)  and  multiple  ranges  (250m  to  2km  while  recent  efforts  have  considered 
ranges  up  to  5km). 

The  “true”  environments  for  these  simulations  (based  on  SW06  test  scenarios)  are  shown  in  Fig.l 
and  Table  1  allowing  for  thin,  medium,  and  thick  sediment  layers.  Consideration  of  a  variety  of 
scenarios  helps  to  address  concerns  that  our  conclusions,  particularly  with  regard  to  frequency,  are 
very  dependent  on  sediment  thickness.  For  the  “exact”  inversion  processing  to  be  discussed  below 
we  shall  assume  that: 

•  the  bottom  consists  of  a  single  linear  sediment  layer  (specified  by  ctop,  7,  and  hsed,  over  a 
half-space  with  sound-speed  chsp )  (parameters  will  vary  depending  on  the  layer  thickness 
under  consideration;  see  Table  1); 

•  all  water  depths  D  will  be  within  78  to  86m  (parameter  value  will  vary  depending  on  the 
layer  thickness  under  consideration;  see  Table  1); 

•  the  fixed  source  ranges  rge  will  each  be  less  than  about  2km; 

•  the  ocean  sound-speed  c(z)  will  vary  with  depth  only  (no  inversion  on  c(z);  see  the  solid 
curve  in  Fig.  1); 

•  Zsou  wiH  he  fixed  (no  inversion  done  on  zsou)\ 

•  we  have  only  one  array  which  will  be  vertical  (VLA)  with  length  56.25m  consisting  of 
16  phones  spaced  at  3.75m  apart  with  array  element  localization  and  top  phone  depth  at 
zph  1  =  14.6m  or  zphi  =  15.6m.  Alternately,  we  have  also  considered  horizontal  arrays 
(HLAs)  with  various  lengths  and  a  variety  of  phones  numbers,  depths  zph,  and  spacings  The 
“true”  arrays  have  no  tilt. 

As  in  the  earlier  simulation  work,  we  will  generate  the  “true”  field  using  the  single  depth-variable 
ocean  sound-speed  profile  seen  as  the  solid  curve  of  Fig.  1  (top),  and  as  before  we  shall  continue 
to  generate  the  synthetic  acoustic  fields  via  RAMGEO  (Collins,  ’94). 

Broadband  (BB)  frequency  averaging  for  the  improvement  of  matched  field  processing  (MFP)  has 
been  around  for  over  twenty  years,  particularly  for  the  suppression  of  sidelobes  (Tolstoy,  ’93). 
For  each  f  the  MFP  values  at  sidelobes  (non-true  parameter  values)  can  vary  quite  a  bit.  Simple 
incoherent  summation  'Pun,ave( p)  (see  Tolstoy,  ’12b)  can  be  quite  stable  with  values  following 
the  high  level  MFP  values  while  not  diminishing  much  for  an  occasional  low  MFP  value.  Such  an 
approach  (with  the  linear  processor)  is  often  used  for  data  plus  a  Bayesian  inversion  to  estimate 
geoacoustic  parameters  (Chapman  &  Jiang,  ’08).  Incoherent  linear  summation  is  also  an  ingredient 
in  many  inversions  using  genetic  algorithms  (Gerstoft  et  al.,  ’03). 
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Figure  1:  Plot  of  simulated  SW06  environment  where  the  upper  subplot  shows  the  ocean  sound- 
speed  c(z )  used  in  the  simulations  (the  exact  and  approximate  profiles  shown),  and  the  lower 
subplot  shows  the  ocean  waveguide  assuming  the  linear  sound-speed  profile  in  a  single  sediment 
layer  over  a  half-space  basement.  Actual  bottom  parameter  values  to  be  found  in  Table  1.  For  cdl 
exact  scenarios  we  will  have  true  zphi  =  14.6m  and  true  ranges  rge  =  265,  480,  780,  980,  and 
2075m. 
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Table  1:  Table  of  geometric  and  environmental  values  for  the  three  sediment  thicknesses  (simu¬ 
lated)  considered.  hse(i  is  the  sediment  thickness,  csed  is  the  sound-speed  at  the  top  of  the  sedi¬ 
ment,  7  is  the  sound-speed  gradient  in  the  sediment  (sound-speed  at  the  bottom  of  the  sediment  is 
given  by  csed  +  7 hsed),  Chsp  is  the  sound-speed  of  the  basement  half-space,  zsou  is  the  source 
depth,  and  D  is  the  water  depth. 


Unfortunately,  this  summation  approach  does  not  indicate  when  component  frequencies  have  con¬ 
tributed  low  levels.  Moreover,  the  summation  MFP  level  can  remain  high  even  when  a  small  subset 
of  components  values  is  legimately  very  low.  Thus,  the  summation  method  can  smear  out  MFP 
levels  so  that  sensitivity  is  actually  reduced  (the  curve  is  less  peaked  with  slower  sidelobe  level 
degradation)  making  it  harder  to  find  “true”  parameter  values.  Thus,  one  can  trade  “robustness” 
for  “sensitivity”.  See  Tolstoy,  ’12b,  for  greater  detail  on  the  method. 

Consider  the  “thin”  sediment  case.  In  Fig.  2  we  see  the  linear  MFP  behavior  as  a  function  of  <7 lsp 
for  a  few  frequencies  (as  indicated:  25,  50,  75,  and  100Hz)  and  at  rge  =  980m.  We  assume  here  that 
all  parameters  other  than  Chsp  are  known  exactly.  First,  we  note  that  there  are  sidelobe  differences 
per  frequency  (as  expected).  Next,  we  note  that  there  are  sensitivity  differences  with  maximum 
sensitivity  at  25Hz.  Finally,  we  note  that  although  the  behavior  varies  as  a  function  of  frequency,  it 
is  impossible  to  predict  the  specific  behavior  as  it  does  not  vary  systematically.  For  our  example, 
frequency  100Hz  shows  more  sensitivity  than  50Hz  while  75Hz  is  more  sensitive  than  100Hz. 


SW06 


thin  sediment  —  exact  @  980m 


chsp  (m/s) 


Figure  2:  Multi-frequency  linear  processor  results  at  a  thin  sediment  for  cj lsp  at  rge  -  980m. 


Consider  a  broadband  average  (over  16  frequencies  25  to  100Hz,  A /  =  5Hz).  In  Fig.  3a  we 
see  that  the  average  'Plin,  ave( p)  (the  solid  line  with  the  filled  in  black  circles)  has  improved  the 
sensitivity  and  reduced  sidelobes  compared  to  some  but  not  all  /,  e.g.,  compared  to  50Hz  but  not 
compared  to  25  or  75Hz  (as  seen  in  Fig.  2).  If  we  consider  Vrninip)  as  shown  by  the  dashed 
line  with  the  open  circles  (and  in  this  example  the  curve  actually  corresponds  to  that  seen  in  Fig. 
2  for  25Hz),  then  we  now  have  significantly  more  sensitivity  than  for  the  straight  average.  That 
is,  we  see  significant  sidelobe  reduction  for  the  new  processor.  Thus,  for  this  parameter  c })sp  we 
have  been  able  to  improve  sensitivity  relative  to  the  averaging  processor,  and  the  new  processor 
is  dominated  automatically  by  the  low  frequency  25Hz  component.  Consider  another  parameter: 
hsed- 

In  Fig.  3b  we  see  the  behavior  of  Pun,ave(p)  versus  Vmin( p)  for  the  more  problematic  parameter 
hsed  where  p  varies  from  (1806,5,1622,2)  to  (1806,55,1622,2)  in  lm  increments  for  hsed .  This 
parameter  often  has  a  number  of  local  maxima  (as  seen  here)  thereby  complicating  convergence 
for  inversion  methods.  First,  we  again  see  that  while  the  averaged  linear  method  shows  some 
sensitivity  to  this  parameter,  this  sensitivity  is  significantly  increased  for  the  minimization  method. 
That  is,  sidelobes  overall  have  been  reduced  much  more  efficiently  for  /Pmin.  Second,  for  this 
parameter  hsed  the  new  processor  has  been  dominated  by  the  85-90Hz  components  (not  shown 
individually  and  we  recall  that  here  c(z)  is  known  exactly).  That  is,  we  again  have  that  the  new 
processor  is  dominated  automatically  by  certain  frequencies. 

In  Fig.  3c  we  see  the  behavior  of  csed  where  the  higher  frequencies  contribute  the  best  resolution 
(80Hz  and  above  where  we  still  have  c(z)  known  exactly).  In  general,  it  seems  that  csed  has  the 
best  sensitivity  at  the  highest  frequencies.  Unfortunately,  those  higher  frequencies  are  most  suscep¬ 
tible  to  experimental  errors.  We  again  notice  that  (as  for  the  other  parameters)  the  minimization 
processor  shows  the  better  sidelobe  reduction  for  our  unknown  parameters,  and  it  is  dominated 
automatically  by  certain  frequencies. 

These  conclusions  still  hold  true  at  longer  ranges,  and  for  other  component  processors  (such  as  the 
minimium  variance  rather  than  the  linear).  That  is,  sensitivity  is  improved  with  the  new  processor 
with  automatic  emphasis  on  the  most  sensitive  /.  Clearly,  the  new  processor  shows  promise  rela¬ 
tive  to  the  averaged  BB  linear  processor  -  when  things  are  known  perfectly.  However,  one  of  the 
strong  points  for  the  averaging  MFP  is  that  when  things  are  not  known  exactly,  i.e.,  when  there  are 
errors  in  our  assumptions  (as  in  a  test  situation),  the  usual  appproach  is  known  to  be  quite  stable 
and  not  overly  disturbed  by  an  occasional  frequency  misstep.  What  happens  for  Vmin ? 

In  particular,  let  us  assume  that: 

•  7  =  0  (thus  we  will  be  inverting  for  the  approximate  cave  in  the  sediment  layer; 

•  source  range  rge  is  known  only  to  within  100m  (rge  and  water  depth  D  are  known  to  be 
linearly  related  -  even  broadband  processing  cannot  separate  out  the  “true”  values  of  rge 
and  D:  see  Tolstoy  et  al.,  ’  12a).  Thus,  we  will  invert  only  for  D  assuming  a  “known”,  i.e., 
approximate,  rge; 

•  the  ocean  sound-speed  c(z)  will  be  assumed  by  the  dashed  (incorrect)  curve  of  Fig.  1; 

•  zsou  wiH  be  fixed  at  30m  (rather  than  its  “true”  value  which  will  vary  as  in  Table  1); 


(a) 


(b) 


SW06 


SW06 


thin  sediment  —  exact  @  980m 


thin  sediment  —  exact  @  980m 


chsp  (m/s) 


SW06 


thin  sediment  —  exact  @  980m 


(c) 


Figure  3:  'Pave.Un  and  'Pmin  (linear  MFP  components)  for  16  frequencies  25  to  100Hz  at  5  Hz 
increments,  with  rge  =  980m,  assuming  the  thin  sediment  scenario  and  other  parameters  known 
exactly,  (a)  The  parameter  considered  is  Cf lsp  varying  from  1700m/s  to  2 100m/ s.  The  true  value  is 
1806m/s.  (b)  The  parameter  considered  is  hsed  varying  from  5  m  to  55  m.  The  true  value  is  12m.  (c) 
The  parameter  considered  is  csed  varying  from  1550m/s  to  1700m/s.  The  true  value  is  1622m/s. 


•  for  the  VLA  we  will  assume  top  phone  depth  at  zph\  =  15.6m  (a  shift  error  of  lm  depth  for 
all  phones). 

In  Fig.  4a  (comparable  to  Fig.  3a  but  at  rge  =  2075m  plus  our  errors)  we  see  that  Vmin( p)  has  de¬ 
graded  a  lot  with  lower  peaks  of  0.71  atp  =  (1800, 12, 1630,0),  0.77  atp  =  (1800, 11, 1630,0), 
and  0.74  at  p  =  (1800, 12, 1640,  0).  We  also  observe  that  non-peak  values  have  again  decreased 
much  more  than  T’ave,iin(p)  ■  and  that  the  new  peaks  can  be  in  different  erroneous  locations  (not 
equivalent  to  simply  raising  the  linear  processor  to  some  power).  Thus,  even  in  the  presence  of 
our  errors  Vmin{ p)  shows  very  good  sidelobe  reduction.  But  what  about  those  rather  low  (and 
incorrectly  shifted)  peak  values?  These  low  values  (less  than  0.9)  turn  out  to  occur  consistently  at 
frequencies  above  60  Hz.  Such  higher  frequency  values  are  understandably  more  affected  by  our 
errors.  However,  if  we  restrict  our  new  processor  to  /  25-60  Hz,  then  we  get  the  results  seen  in 
Fig.  4  for  the  starred  curves,  i.e.,  great  performance  with  correct,  strong  peaks  and  much  reduced 
sidelobe  s. 

Why  do  we  care  about  this  proposed  method?  First,  the  method  offers  improved  efficiency  in 
sequential  inversion  computations.  In  particular,  let  us  begin  the  inversions  at  a  parameter  point  p 
at  a  low  frequency  (sequential  with  frequency)  and  then  increase  frequency  in  steps.  If  at  any  step 
we  find  that  'Plin  (/,  p)  is  less  than  our  target  amount  (a  fairly  high  value  set  by  the  user),  then  we 
can  cease  the  frequency  computations  and  move  to  another  parameter  point.  Thus,  many  points 
may  be  quickly  eliminated. 

Additionally,  at  LFs  we  can  employ  cruder  sampling  of  the  multi-dimensional  solution  space.  Re¬ 
stricting  the  processor  to  rather  low  frequencies  such  as  25-60Hz  can  mean  that  fewer  parameters 
at  cruder  sampling  intervals  will  need  to  considered  to  find  all  the  peaks. 

This  method  also  indicates  when  frequency  difficulties  appear.  That  is,  consistently  low  MFP 
values  suggest  a  problem.  Then,  the  user  can  study  why  and  when  such  persistent  low  MFP  values 
occur.  Is  the  signal-to-noise  level  low?  Is  there  an  erroneous  asssumption  somewhere?  In  general, 
such  results  can  act  as  a  helpful  source  of  debugging  as  well  as  of  general  information  about  the 
inversion  itself. 

Next,  this  new  method  can  be  insensitive  to  expected  errors  in  parameters  assumed  to  be  “fixed” 
such  as  zphi .  The  technique  allows  for  more  error  in  the  experimental  measurements. 

Finally,  this  method  also  offers  improved  parameter  resolution.  In  particular,  it  automatically 
emphasizes  those  frequencies  which  are  sensitive  to  our  parameters  of  interest.  That  is,  it  does 
not  smear  out  that  sensitivity  over  a  wide  range  of  frequencies  (as  happens  with  general  incoherent 
summation)  but  rather  is  dominated  by  the  sensitivity  offered  by  any  frequency  component. 

What  are  potential  problems?  This  approach  will  not  work  if  there  are  low  signal-to-noise  levels 
at  any  of  the  final  component  frequencies.  A  requirement  for  this  method  is  that  each  frequency 
contribution  not  be  “too”  erroneous.  A  band  of  frequencies,  e.g.,  above  60Hz,  may  be  eliminated  if 
their  contributions  are  problematic  -  so  long  as  those  difficulties  are  understood  (severely  affected 
by  expected  experimental  errors  or  false  assumptions  about  the  environment?). 

Thus,  we  conclude  from  this  year’s  work  that  the  new  minimization  method: 


•  is  broadband; 
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Figure  4:  Same  as  Fig.  3  except  that  small  errors  have  been  assumed  (environmental  parameters 
are  known  only  approximately )  and  we  are  at  rge  =  2075m.  We  also  have  curves  (starred,  dashed) 
for  'Prni  n  computed  using  only  frequencies  25-60Hz  (we  recall  that  the  low  peak  values  ended  up 
indicating  errors  in  the  higher  frequencies),  (a)  Chsp,  (b)  hsed,  and  (c)  csed. 


•  will  be  dominated  automatically  by  those  frequencies  which  are  most  sensitive  to  unknown 
parameters; 

•  improves  sensitivity  and  resolution  compared  to  incoherent  broadband  summation  (for  any 
component  processor  such  as  the  linear  or  minimum  variance  processor); 

•  can  improve  inversion  efficiency  (for  sequential  /  computations); 

•  shows  sidelobe  improvement  even  in  the  presence  of  (our  expected)  errors; 

•  indicates  when  something  is  wrong,  i.e.,  the  peak  value  will  be  “low”.  Such  problems  need 
to  be  pursued  to  understand  why  that  has  happened.  The  method  can  also  indicate  when  a 
band  of  component  frequencies  has  a  problem,  and  this  can  lead  to  the  adoption  of  only  the 
low  frequency  components. 


As  a  final  note,  implementation  of  the  new  method  should  be  easy  -  it  can  simply  replace  the 
incoherent  summation  of  components  with  a  minimization  function. 

WORK  COMPLETED 

Recent  work  (FY08)  completed  includes: 

•  continued  development  of  a  new  BB  signal  processing  method  0Pmin(/)); 

•  application  of  the  minimization  method  to  several  simulated  SW06  scenarios  (single  sedi¬ 
ment  layer:  thin,  medium,  or  thick); 

•  extension  to  longer  ranges  (incorporation  of  KRAKEN  for  more  speed  assuming  range- 
independence)  in  anticipation  of  application  to  Knobles  data; 

•  examination  of  the  minimization  method  sensitivity  as  a  function  expected  errors  (in  ocean 
sound-speed  c(z),  source  depth  zsou,  and  in  VLA  depth  zphi)- 

RESULTS 

We  have  a  new  BB  processor  'Pmin(f)  which  promises  excellent  resolution  for  G.I.  at  LFs  (/ 
within  25-75Hz)  and  at  close  ranges  (within  250m  to  about  1km),  and  even  in  the  presence  of 
expected  test  errors.  This  processor  will  also  indicate  when  it  has  trouble  (it  will  show  low  MFP 
values). 

IMPACT/APPLICATION 

As  a  result  of  the  work  this  past  year  we  have  developed  and  better  understand: 

•  the  FF  G.I.  method  as  applied  to  a  variety  of  simulated  SW06  data,  particularly  with  regard 
to  sensitivity  for  bottom  parameters  as  a  function  of  frequency  and  range; 

•  the  effects  of  expected  errors  in  a  SW06  test  environment; 

•  a  new  BB  inversion  method  (relative  to  standard  BB  incoherent  averaging)  with  demon¬ 
strated  success  on  simulated  SW06  data; 


•  the  potential  success  of  HLAs  for  G.I. 


RELATED  PROJECTS 

The  G.I.  work  is  related  to  work  by  R.  Chapman  and  colleagues  (U.  Victoria),  D.  Knobles  and 
colleagues  (U.  Texas  at  Austin),  W.  Hodgkiss  and  colleagues  (Scripps),  and  other  researchers  in 
SW06  and  shallow  water  inversion  (such  as  P.  Gerstoft,  P.  Nielsen,  C.  Harrison). 
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